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Abstract 

We use ab initio calculations to systematically study the phase diagram of multiferroic 
Sri_ 3 ;Ba 3 ;Mn 03 (0 < x < 1) as a function of chemical doping, epitaxial strain and external pres¬ 
sure. We find that by replacing Sr with Ba in cubic SrMnOs and imposing epitaxial strain, the 
material can be tuned to the vicinity of a first order transition between two multiferroic phases, one 
antiferromagnetic with a smaller polarization and one ferromagnetic with a larger polarization. A 
giant effective magneto-electric coupling and cross-field control (electric field control of magnetism 
or magnetic field control of polarization) can be achieved in the vicinity of the transition. The 
dependence of the theoretically computed transition point on the choice of exchange correlation 
functionals is determined and is found to be non-negligible. We also show that the perovskite 
structure of BaMnOs can be stabilized relative to its hexagonal polymorphs at pressures larger 
than 20 GPa. 
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I. INTRODUCTION 


Searching for multiferroic materials, in particular those simultaneously possessing fer¬ 
roelectric polarization and ferromagnetic moments with strong coupling between the two 
order parameters, is currently of great interest [l|, 2 |. Among all single-phase multiferroic 
materials, perovskite Sri_a,B^Mn03 (0 < x < 1 ) stands out as a promising candidate 


for multiferroics in theory 


8 |. Ferroelectricity was theoretically predicted in perovskite 


BaMnOs j^. It was also predicted [w, 11 | that under biaxial strain, perovskite SrMnOs (a 
paraelectric antiferromagnet at ambient conditions) becomes a ferroelectric ferromagnet. 

In this paper, we examine the phase diagram of Sri_3;Ba3;Mn03 (0 < x < 1 ) as a func¬ 
tion of chemical doping, epitaxial strain and external pressure. We show that by doping 
cubic SrMn03 with an appropriate concentration of Ba and imposing epitaxial strain, the 
material can be tuned to the vicinity of a hrst order phase transition between two multifer¬ 
roic phases, one antiferromagnetic with a smaller polarization and one ferromagnetic with 
a larger polarization. In the vicinity of the transition point, reasonable electric helds can 
be used to switch the material between ferromagnetic and antiferromagnetic states while 
feasible magnetic helds can change the electrical polarization. However, the critical Ba dop¬ 
ing and the critical epitaxial strain strongly depend on the exchange correlation functionals. 


Building on our previous work 


12 | . we compare three different versions of Perdew-Burke- 


Ernzerhof (PBE) parameterized exchange correlation functionals and show that they make 
substantially different predictions for the position of the magnetic transition boundary in 
the phase diagram of Ba doping and epitaxial strain. We focus on the perovskite structure of 
Sri_2.Ba3;Mn03 to achieve effective magneto-electric couplings. However, BaMn03 can also 
crystallize in hexagonal polymorph, which is more stable than the perovskite structure. We 
show that applying an external pressure of over 20 GPa on BaMnOs is a feasible approach 
to stabilize the more interesting perovskite structure over its hexagonal polymorph. 

The rest of this paper is organized as follows. Sec. UTI presents the computational details. 
Sec.im discusses the phase diagram of perovskite Sri-ajBa^MnOs as a function of chemical 
doping and estimates the helds required for magneto-electric cross-held control. Sec. HVl 
presents the magnetic phase diagram of perovskite SrMn03, Sro.5Bao,5Mn03 and BaMn03 
as a function of strain. Sec. El studies energetic dependence of perovskite and hexagonal 
polymorphs of BaMn03 on external pressure. We conclude in Sec. IVII 
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II. COMPUTATIONAL DETAILS 


We perform density functional theory calculations 


wave approach 


3, Cl 


within the ab initio plane- 


15| . as implemented in the Vienna Ab-initio Simulation Package (VASP) 16 |. 
We employ projector augmented wave (PAW) pseudopotentials [l7, 18|. We use an energy 
cutoff 600 eV. All the calculations allow for spin-polarization to study different types of 
long-range magnetic orderings. For most calculations, both cell and internal coordinates 
are fully relaxed until each force component is smaller than 10 meV/A and stress tensor 
is smaller than 10 kBar. For the study of strain, the in-plane lattice constants are fixed 
during atomic relaxation; the internal coordinates and out-of-plane lattice constant are fully 
relaxed. For the study of pressure, the atomic relaxation is converged to a specihc external 
pressure. 

We consider different structures of Sri_ 2 .Ba,j,Mn 03 . A 12 x 12 x 10 Monkhorst-Pack grid is 
used to sample the Brillouin zone of the simulation cell used in Sec. Illll and Sec. llVI to study 
the perovskite structure. For the hexagonal polymorphs studied in Sec. IVl a 10 x 10 x 10 
Monkhorst-Pack grid is used. The simulation cells are provided in each section for ease of 
reading. 

Our recent study shows |3] that different exchange correlation functionals make sub¬ 
stantially different predictions for structural and magnetic properties of ferroelectric man- 
ganites. By comparing to available experimental data for Sro.sBao.sMnOs alloy, we establish 
that charge-density-only generalized gradient approximation with Perdew-Burke-Ernzerhof 


parameterization 


lOj plus Hubbard U and Hund’s J corrections (PBE-|-t/-|-J) |20j and spin- 


polarized generalized gradient a ppr oximation with Perdew-Burke-Ernzerhof parameteriza¬ 
tion revised for solids (sPBEsol) 2l|] make overall better predictions than the spin-polarized 
PBE (sPBE) and local spin density approximation (LSDA). In this paper, for most of the 
results, we use PBE+U+J and sPBEsol. For reference, we also compare PBE-I-17-I-J and 
sPBEsol to sPBE plus effective Hubbard Ues corrections. For the PBE-I-17-I-J method, we 
use accepted values of U = 5 eV and J = 0.7 eV for Mn d orbitals unless otherwise specihed. 
This method is implemented in VASP by turning on LDAUTYPE = 4. For the sPBE-|-[/eff 
method, we follow previous work |3] and use Ues = 1.7 eV. This method is implemented in 
VASP by turning on LDAUTYPE = 1. 
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FIG. 1: Atomic structure of perovskite BaMnOs. The orange, purple and red balls are Ba, Mn 
and O atoms, respectively. The Mn off-center displacement (^Mn-O is formally defined as: ^Mn-O = 
B -^Mn-o where Run-O are the vectors connecting a Mn ion to its six nearest oxygen neighbors. 
A) Mn off-center displacement 6 along the [100] direction. B) Mn off-center displacement 6 along 
the [110] direction. 


III. CHEMICAL DOPING 


Following our previous work 


12|, we first use the PBE+U+J method to examine the Ba- 


end member compound perovskite BaMnOa, considering four different magnetic orderings 
that are commonly found in manganites: ferromagnetic (F), two-sublattice Neel ordering 
(G-type), (1, 0, 0) stripe ordering (A-type) and (1, 1, 0) stripe ordering (G-type). Fig. [1] 
shows the atomic structure of perovskite BaMnOs. We use a \/2 x \/2 x 2 simulation cell that 
can accommodate all the above four magnetic orderings. In Fig. [H 6 is the Mn off-center 
displacement which is along the [100] axis (Fig. [T]A) and along the [110] axis (Fig. [T]B). 

Fig.[2]presents the dependence of the total energy on the Mn off-center displacement (5Mn-o 
for both [100] and [110] directions. The reference energy is that of the G-type antiferromag- 
net at (5Mn-o = 0. All the other magnetic states exhibit an energy minimum at (Imii-o 7^ 0, 
but the ferromagnetic state is lowest for both directions, and has the largest f^Mn-o at the 
minimum. Interestingly, the energy difference between different magnetic phases depends on 
the direction of polarization. For (5Mn-o along the [110] direction, ferromagnetism is clearly 
more stable than any of the antiferromagnetic states, while for ^Mn-o along the [100] di- 
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FIG. 2: Change in total energy of perovskite BaMnOs as a function of Mn off-center displacement 
(^Mn-O computed using the PBE-|-C/-|-J method with U = 5 eV and J = 0.7 eV for ferromagnetic 
{F, red curves), A-type antiferromagnetic [A, orange curves), C-type antiferromagnetic (C, green 
curves) and G-type antiferromagnetic (G, blue curves) states (for the definition of magnetic or¬ 
derings, see the main text) for A) along the [100] and B) along the [110] directions of Jmii-o- 
The zero energy state is chosen as the energy of the G-type antiferromagnetic state at ^Mn-O = 0. 
Calculations were hrst performed using a cubic structure (^Mn-O = 0) and a fully relaxed structure, 
then were interpolated by using intermediate structures. 


rection, ferromagnetism and A-type antiferromagnetism are extremely close in energy. The 
near coincidence of energies suggests that for x slightly less than 1, it may be possible to 
switch the magnetic states by applying an electric field and rotating the orientation of ^Mn-o- 
The results shown in Fig. [3] establish that within the PBE-I-17-I-J method with U = 5 
eV and J = 0.7 eV, perovskite BaMnOa is a ferromagnet with strong Mn off-center dis¬ 
placements. Now we consider the evolution of magnetic and ferroelectric states with Ba 
concentration x in Sri_ 3 .Baa;Mn 03 . Fig. [3] presents the dependence of the total energy of 
two magnetic states of Sri_ 3 ;Baa;Mn 03 (for x = 0.0, 0.25, 0.5, 0.75 and 1.0) on the inversion- 
symmetry-breaking Mn-0 displacement i^Mn-o along [100] and [110] directions. The hgure 
shows that for both directions of ^Mn-o, the G-type antiferromagnetic state (which is insu¬ 
lating for all x) has a paraelectric-to-ferroelectric transition as the Ba doping x increases, 
consistent with previous studies jd-OjIo]. We also see that for all x the ferromagnetic state 
always has a nonzero Mn off-center displacement ^Mn-o which is larger than that of antifer¬ 
romagnetic state at the same x. With increasing Ba doping the energy difference between 
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FIG. 3: Total energy of Sri_ 3 ;Baa;Mn 03 (0 < x < 1) calculated using the PBE+f7+J method 
with [7 = 5 eV and J = 0.7 eV as a function of Mn off-center displacement i^Mn-o- The solid lines 
indicate ferromagnetic ordering and the dashed lines indicate G-type antiferromagnetic ordering. 
The red, orange, green, blue and indigo curves correspond to Ba doping of x = 0, 0.25, 0.5, 0.75 and 
1, respectively. For each x, the zero of energy is chosen as the energy of the G-type antiferromagnet 
at (^Mn-o = 0. A) ^Mn-O along the [100] direction; B) (5Mn-o along the [110] direction. 

the ferromagnetic and antiferromagnetic states rapidly decreases and for the parameters U 
= 5 eV and J = 0.7 eV used here, a hrst-order antiferromagnetic-to-ferromagnetic tran¬ 
sition occurs around x = 0.75 (blue curves), accompanied by an abrupt jump in the Mn 
off-center displacement. Therefore, in the vicinity of the transition, we can use an external 
electric held to increase the Mn off-center displacement and change the magnetic ordering, 
or use an external magnetic held to induce a hnite change in Mn oh-center displacements. 
This ehective magneto-electric cross-held control is similar to what is found in epitaxially 
strained SrMnOs, indicating that chemical doping is a complementary approach to tuning 
multiferroic properties. 

Based on Fig. [3l we can make some estimation of critical helds required for switching. 
For physically reasonable parameters U = 5 eV and J = 0.7 eV, our calculations show 
that the critical doping x is around 0.75. At x = 0.75 and along the [110] direction, the 
two lowest energy states are G-type antiferromagnetism (G) and ferromagnetism (F). The 
critical helds are estimated by: 
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FIG. 4: Energy difference between ferromagnetic state {F) and G-type antiferromagnetic state 
(G) of Sri_ 3 ;Baa;Mn 03 (0 < x < 1) calculated using difference exchange correlation functionals: 
charge-density-only PBE functional with Hubbard U and Hund’s J correction (PBE-|-17-|-J) with 
U = 5 eV and J = 0.7 eV; spin polarized PBE for solids (sPBEsol) and spin polarized PBE with 
effective Hubbard Uqs correction (sPBE-|-17efr). A) ^Mn-O along the [100] direction; B) Jmii-o along 
the [110] direction. 


■ {Pf - Pg) ^ E{F) - E{G) 


( 1 ) 


7^, • M ~ E{F) - E{G) (2) 

where Pp = 56 /iC/cm^ and Pq = 26 /rC/cm^ are electric polarizations calculated for 
ferromagnetic and G-type antiferromagnetic states with on-site magnetic moment M ~ 3/i^ 
using the Berry phase method |22t] (their electric held dependence is neglected). We hnd 
that £c is about 42 MV/m and Fie is about 25 T. While these helds are large, they are 
experimentally achievable and have been applied to, for example, BaTiOs thin hlms 23|. 
More importantly, tuning the Ba doping x by a few percent can rapidly reduce these critical 
helds. Use of other interaction parameters U = A eV and J = 0.6 eV (the smallest U value 
to stabilize a high-spin state in cubic SrMnOa using the PBE-|-[/-|-J method 2^) shows that 
at X = 1, we can still switch the two lowest magnetic states by a critical magnetic held of 
12 T or switch those two polarizations by an electric held of 30 MV/m. As we see, properly 
adjusting the Ba doping x can compensate for the uncertainty of U and J parameters in 
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theory. 

Our previous works [^, 241 show that magnetic properties of transition metal oxides 
depend on the choice of exchange-correlation functionals. To test whether our prediction of 
a doping-driven magnetic transition in Sri_j;Ba 3 ;Mn 03 is robust to the choice of exchange- 
correlation potential, we compare three different types of PBE exchange correlation func¬ 
tionals: spin-polarized PBE revised for solids (sPBEsol), charge-density-only PBE plus Hub¬ 
bard U and Hund’s J correction (PBE-|-f/-|-J) and spin-polarized PBE plus effective f4fr 
correction (sPBE-|-f/efr) by calculating the energy difference of Sri_ 3 ,Ba 3 ;Mn 03 between the 
ferromagnetic and G-type antiferromagnetic states. The comparison is shown in Fig. |H We 
hnd that while PBE+U+J and sPBE-|-f/eff predict that there is an antiferromagnetic to fer¬ 
romagnetic transition as x is increased from 0 to 1 for perovskite Sri_ 3 ;Baa;Mn 03 , sPBEsol 
predicts that perovskite Sri_a;Ba 2 ;Mn 03 is G-type antiferromagnetic for the whole range of 
X (0 < x < 1). 

We can test the functionals by comparing to experimental data for Sro. 5 Bao. 5 Mn 03 al¬ 
loy jj]. We find from our previous study that both PBE-I-17-I-J and sPBEsol make good 
predictions of structural and magnetic properties (better than sPBE and LSDA). Hence 
we can not definitely predict the ground state magnetic ordering of perovskite BaMnOs. 
However, that doping cubic SrMn 03 with Ba favors ferromagnetism is a robust result. If 
the end member perovskite BaMn 03 were ferromagnetic, then the antiferromagnetic-to- 
ferromagnetic transition would occur at an intermediate doping Xc, in the vicinity of which 
an effective giant magneto-electric coupling can be induced. If perovskite BaMn 03 were not 
ferromagnetic, then epitaxial strain may be used to stabilize the ferromagnetic state. 


IV. EPITAXIAL STRAIN 

Strain engineering has been widely used to induce multiferroic states [nl, . Lee and 
Rabe showed that epitaxial strain can induce a ferroelectric-ferromagnetic ground state 
in an otherwise antiferromagnetic cubic SrMnOs with a critical tensile strain of 3.4% (pre¬ 
dicted by sPBE plus an effective Hubbard Ues extension, for short sPBE-|-f/eff). In this 
section, we define the strain A as: 


A 


a — ttQ 
Oo 
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X 100% 


(3) 









where a is a common lattice constant that is imposed on all the magnetic orderings and ao 
is the theoretically calculated lattice constant of the G-type antiferromagnetic state. We 
note that because we define strain with respect to a particular state (the G-type antiferro- 
magnet), the critical strain depends both on the energy difference between ferromagnetic 
and antiferromagnetic states at equilibrium (denoted by 6E) and on the difference between 
equilibrium lattice constant (denoted by 6X). We first compare the predictions of sPBE+f/eg 
to PBE+[/+J and sPBEsol methods for the phase diagram of strained SrMnOa. We show 
that the three methods make substantially different predictions on the critical strain that 
stabilizes the ferroelectric-ferromagnetic ground state of SrMnOs. 

Next we study the phase diagram of epitaxially strained perovskite Sro.sBao.sMnOa alloy 
and epitaxially strained perovskite BaMnOs. We show that Ba doping acts like an “effective 
strain”, so that increasing the Ba concentration in Sri_a;Baa;Mn 03 can reduce the critical 
strain that stabilizes the ferromagnetic state. However, just like epitaxially strained SrMnOs, 
the critical strain as a function of Ba doping also strongly depends on exchange correlation 
functionals. 


Perovskite SrMnOq 


Fig. Elshows the total energies of SrMnOs and Mn off-center displacement along the [100] 
direction under tensile strain. Fig. [5]A1 and B1 are calculated using the sPBEsol method. 
Fig. [5]A.2 and B2 are calculated using the PBE-I-G-I-J method with U = 5 eV and J = 0.7 
eV. Fig. [5]A.3 and B3 are calculated using the sPBE-|-17eff method with f/gff = U—J = 1.7 eV, 
following Ref. lOj. The panels A of Fig. [5] show the total energy of tensile strained SrMnOs 
with different magnetic orderings (F, A, G and G, as dehned in the previous section). 
We find the predicted critical tensile strain needed to drive a ferromagnetic transition is 
6.8% for sPBEsol (not explicitly shown in the Fig. [5]A1), 4.8% for PBE-|-f/-|-J and 3.2% for 
sPBE-|-7/pff , respectively. The critical strain for sPBE-|-f/eff is slightly smaller than that given 


in Ref. 


lOj (see 2^). We notice that all the three methods predict that with increasing 


tensile strain, the elastic energy eventually stabilizes the ferromagnetic state over the G- 
type antiferromagnetic state. However, the strain dependence are substantially different. 
The difference has two origins. One is the energy difference 5E between the ferromagnetic 
and G-type antiferromagnetic states at equilibrium, which is shown as the dashed vertical 
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SrMnOs 



FIG. 5: Perovskite SrMnOs; Al) and Bl) using the sPBEsol method. Al) Total energy for 
different magnetic states under tensile strain. The zero point is chosen as the total energy of 
G-type antiferromagnetic state at zero strain. A is the strain with respect to the equilibrium 
lattice constant of G-type antiferromagnetic state. Bl) Mn off-center displacement along the [110] 
direction of SrMnOs for different magnetic states under tensile strain. A2) and B2) using the 
PBE-|-?7-|-J method with U = 5 eV and J = 0.7 eV. A2) same as Al), B2) same as Bl). A3) and 
B3) using the sPBE-l-Gefr method with Geff = U — J = 1.7 eV. A3) same as Al), B3) same as Bl). 
The red squares and lines are the ferromagnetic state {F). The orange squares and lines are the 
A-type antiferromagnetic state {A). The green squares and lines are the G-type antiferromagnetic 
state (G). The blue squares and lines are the G-type antiferromagnetic state (G). The horizontal 
dashed line denotes dA, which is the difference in equilibrium lattice constants between F and G 
orderings. The vertical dashed line denotes 5E, which is the equilibrium energy difference between 
F and G orderings. 
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FIG. 6: Calculations of total energies of SrMnOs for different magnetic states under tensile strain, 
using the PBE+I/+J method (A1 U = 4 eV; A2 U = 5 eV; A3 U = 6 eV; for all three cases, 
J = 0.7 eV). The red squares and lines are the ferromagnetic state {F). The orange squares 
and lines are the A-type antiferromagnetic state (A). The green squares and lines are the C-type 
antiferromagnetic state (C). The blue squares and lines are the G-type antiferromagnetic state 
(G). The horizontal dashed line denotes 5X, which is the difference in equilibrium lattice constants 
between F and G orderings. The vertical dashed line denotes 5E, which is the equilibrium energy 
difference between F and G orderings. 

lines in the panels A of Fig. |5l The other is the difference in equilibrium lattice constants 
6X between the ferromagnetic state and G-type antiferromagnetic state, which is shown 
as the horizontal lines in the panels A of Fig. |5l 6E is more important. Increasing 6E 
substantially increases the critical strain. sPBEsol, PBE-I-17-I-J and sPBE-t-f/eg predict 6E 
to be 149 meV/Mn, 83 meV/Mn and 57 meV/Mn, respectively. For a hxed 6E, increasing dX 
(within a reasonable range) can reduce the critical strain, because at the same constrained 
lattice constant, a larger 6X means that the ferromagnetic state is under smaller strain, thus 
inducing smaller elastic energy cost. The sPBEsol, PBE+U+J and sPBE-|-[/efr methods 
predict that the equilibrium lattice constant of the ferromagnetic state is 1.4%, 0.5% and 
1.6% larger than that of the G-type antiferromagnetic state. sPBE-t-Geff predicts the smallest 
dE and the largest 6X, the combination of which leads to the smallest critical strain (3.4%). 
On the other hand, sPBEsol predicts the largest 6E, which leads to the largest critical strain 
( 6 . 8 %). 

The panels B of Fig.[5]show the Mn off-center displacement ^Mn-o along the [110] direction 
for strained SrMnOs. We hnd that all the three methods yield similar behavior for (Jmii-o 
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FIG. 7: Atomic structure of perovskite Sro.sBao.sMnOs. The orange, green, purple and red balls 
are Sr, Ba, Mn and O atoms, respectively. A) Mn off-center displacement 6 along the [100] axis 
direction. B) Mn off-center displacement 6 along the [110] axis direction. For tensile strain, we 
focus on Mn off-center displacement 6 along the [110] axis direction. 


with sPBE-|-t4fr predicting a slightly larger displacement than sPBEsol and PBE+U+J, 
consistent with our previous study [12]. 

Finally we briefly discuss the Hubbard U dependence. We study the U dependence using 
the PBE+U+J method, but the conclusion also applies to sPBE-|-Hefr. We study seyeral 
different values of Hubbard U and compare the results in Fig. [6l As we increase the Hubbard 
U from 4 to 5 to 6 eV, 6E (dashed vertical lines in Fig. E]) is reduced from 103 meV/Mn to 
83 meV/Mn to 66 meV/Mn. At the same time, 6X (dashed horizontal lines in Fig. [6]) also 
decreases from 1.2%, to 0.5% to 0.4%. While a decreasing 6E favors the ferromagnetism and 
a decreasing 6X disfavors the ferromagnetism, 6E is the dominating factor and eventually 
the critical strain is reduced. Therefore, increasing U favors the ferromagnetism, which 
is consistent with the simple picture that a larger U suppresses the superexchange and 
antiferromagnetism, and thus ferromagnetism is favored. 


B. Perovskite Sro.sBao.sMnOs alloy 

In this sub-section, we study perovskite Sro.sBao.sMnOs alloy under epitaxial strain. To 
simulate cation alloys, we employ a simulation cell with a checkerboard arrangement of 
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Sr-Ba cations (see Fig. [^. 

To be consistent with the study of strained SrMnOs in the previous sub-section, we 
consider the Mn off-center displacement along the [110] direction and apply tensile strain, 
which can enhance the [110] polarization. Following the same convention, the biaxial strain 
A is defined in Eq. ([3]). Like SrMnOs, we also study four common magnetic orderings {F, 
A, C, G) in Sro.sBao.sMnOs alloy. The panels A of Fig. |8] present the total energy of 
different magnetic orderings of the Sro.sBao.sMnOs alloy as a function of tensile strain, using 
sPBEsol, PBE+U+J and sPBE-t-Fefr- First we notice that for all three methods, at A = 
0, the G-type antiferromagnetic state has lower energy than the ferromagnetic state. With 
an increasing lattice constant a, the total energy of the G-type antiferromagnetic ordering 
monotonically increases, while the total energy of the ferromagnetic state first decreases to 
the minimum at its own equilibrium position and then increases with strain. In particular, 
we notice that for the sPBE-|-Geff method, the minimum of the ferromagnetic state is lower 
than the minimum of the G-type antiferromagnetic state. This means that without any 
epitaxial strain, sPBE-|-f/eff predicts that the ground state of the Sro.sBao.sMnOs alloy is 
ferromagnetic. However, sPBEsol and PBE+U+J predict that the ground state of the 
Sro.sBao.sMnOs alloy is G-type antiferromagnetic. Therefore, using sPBEsol and PBE-l-17-1-J 
methods, we predict that we need a critical tensile strain (5.2% and 3.5%, respectively) to 
stabilize the ferromagnetic state. 

For all three methods, the critical strain is predicted to be smaller than that of pure 
SrMnOs because the presence of Ba increases the Mn off-center displacement and favors 
the ferromagnetic ordering |^, which acts as an “effective strain”. The panels B of Fig. [E] 
show the Mn off-center displacement ^Mn-O; which monotonically increases as a function of 
strain. For each method, at the critical strain, the Mn off-center displacement ^Mn-o of the 
ferromagnetic state exceeds 0.25 A, which is larger than the Ti off-center displacements of 
commercially important ferroelectrics BaTiOs. 

C. Perovskite BaMnOs 

In this sub-section, we study perovskite BaMnOs under epitaxial strain. The simulation 
cell is shown in Fig. [T] We follow the convention of the previous two sub-sections. The 
results are shown in Fig. |9l From the panels A of Fig. [9], both PBE+U+J and sPBE-|-[/eff 
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FIG. 8: Perovskite Sro.sBao.sMnOs: Al) and Bl) using the sPBEsol method. Al) Total energy 
for different magnetic states under tensile strain. The zero point is chosen as the total energy 
of G-type antiferromagnetic state at zero strain. A is the strain with respect to the equilibrium 
lattice constant of G-type antiferromagnetic state. Bl) Mn off-center displacement along the [110] 
direction of SrMnOs for different magnetic states under tensile strain. A2) and B2) using the 
PBE+f7-|-J method with U = 5 eV and J = 0.7 eV. A2) same as Al). B2) same as Bl). A3) 
and B3) using the sPBE+Geff method with t4ff = U — J = 1.7 eV. A3) same as Al). B3) 
same as Bl). The red squares and lines are the ferromagnetic state (F). The orange squares 
and lines are the A-type antiferromagnetic state (A). The green squares and lines are the G-type 
antiferromagnetic state (G). The blue squares and lines are the G-type antiferromagnetic state 
(G). 
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BaMnOs 



FIG. 9: Perovskite BaMnOs: Al) and Bl) using the sPBEsol method. Al) Total energy for 
different magnetic states under tensile strain. The zero point is chosen as the total energy of 
G-type antiferromagnetic state at zero strain. A is the strain with respect to the equilibrium 
lattice constant of G-type antiferromagnetic state. Bl) Mn off-center displacement along the [110] 
direction of SrMnOa for different magnetic states under tensile strain. A2) and B2) using the 
PBE-|-t/+J method with U = 5 eV and J = 0.7 eV. A2) same as Al). B2) same as Bl). A3) 
and B3) using the sPBE-|-t/eff method with Uqs = U — J = 1.7 eV. A3) same as Al). B3) 
same as Bl). The red squares and lines are the ferromagnetic state (F). The orange squares 
and lines are the A-type antiferromagnetic state (A). The green squares and lines are the G-type 
antiferromagnetic state (G). The blue squares and lines are the G-type antiferromagnetic state 
(G). 
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predict that without any epitaxial strain, the ground state of perovskite BaMnOs is ferro¬ 
magnetic. However, sPBEsol predict that the ground state of perovskite BaMnOs is still 
G-type antiferromagnetic ordering and 4% tensile strain is needed to stabilize the ferromag¬ 
netic state. From the panels B of Fig. [HI we hnd that within a 4% tensile strain, all the 
three methods predict a robust Mn off-center displacement for the four magnetic orderings 
investigated. 

Finally, we summarize the critical tensile strain for Sri_a;Baj;Mn 03 (x = 0, 0.5 and 1) 
calculated with different exchange correlation functionals. It is evident that for all three 
methods used here, increasing Ba doping decreases the critical strain that stabilizes the 
ferromagnetic ordering. For each Ba doping x in Sri_ 3 ;Ba 3 ;Mn 03 , sPBE-|-f/eff most favors 
ferromagnetism, while sPBEsol least favors ferromagnetism. This is consistent with Fig. 01 

We also studied the compressive strain (not shown here), which enhances the polarization 
along the [100] direction. For compressive strain, the magnetic energy dependence on strain 
and exchange correlation functional is very similar to what we have found for tensile strain. 


TABLE I: Critical tensile strain that stabilizes ferromagnetic state in Sri_ 3 ;Ba 3 ;Mn 03 (0 < x < 1), 
calculated using the sPBEsol, PBE-|-t/-|-J and sPBE-|-17eff methods. In the PBE-|-17-|-J method, 
we use t/ = 5 eV and J = 0.7 eV. In the sPBE-l-Leff method, we use t/eff = 1-7 eV. 



sPBEsol 

PBE+U+J 

sPBE-|-t/eff 

SrMnOs 

6.8% 

4.8% 

3.2% 

Sro.sBao.sMnOs 

5.2% 

3.5% 

no strain needed 

BaMnOs 

4.0% 

no strain needed 

no strain needed 


V. PRESSURE STABILIZATION OF PEROVSKITE BaMnO. 


While our calculations predict that perovskite Sri_ 3 ;Baa;Mn 03 has appealing multiferroic 
properties, a hexagonal polymorph of manganites has lower energy than the perovskite 
structure 


26 


27l |. However, the perovskite polymorph may be stabilized by pressure [28 1. 


In this section, we focus on BaMn 03 and study the energetics of different polymorphs as a 
function of pressure. All the results presented are obtained using the PBE-I-17-I-J method 
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FIG. 10: Five different atomic structures of BaMnOs. The green balls are La atoms. The pur¬ 
ple cages are MnOg. A) the hexagonal structure with face-sharing MnOe oxygen octahedra; B) 
the hexagonal structure with mixed face-sharing and corner sharing MnOg oxygen octahedra; C) 
the cubic perovskite structure with corner-sharing MnOg oxygen octahedra; D) the tetragonal 
perovskite structure with Mn off-center displacement along the z direction; E) the orthorhombic 
perovskite structure with Mn off-center displacement in the xy plane. The orange arrows indicate 
spin arrangements in the antiferromagnetic state considered for each structure. 

{U = 5 eV and J = 0.7 eV). Other methods, such as sPBEsol and sPBE, yield qualitatively 
consistent conclusions, but with a critical pressure 20% larger in magnitude (see below). 

Fig. [To] shows the atomic structure of hve different polymorphs of BaMnOa. Fig. [TOlA 
shows the structure of hexagonal BaMnOa with face-sharing MnOg oxygen octahedra. This 
structure is called 2H. Fig. [TUB is also a hexagonal structure of BaMnOa with mixed face¬ 
sharing and corner-sharing MnOg oxygen octahedra. This structure is called AH. Fig. [TUjC is 
the cubic perovskite structure of BaMnOa with corner-sharing MnOg octahedra. This struc¬ 
ture is referred to as “cubic”. Fig. HOD is a tetragonal perovskite structure of BaMnOa with 
Mn off-center displacements along the z direction. This structure is referred to as “tetrag¬ 
onal”. Fig. [TUB is an orthorhombic perovskite structure of BaMnOa with Mn off-center 
displacements in the xy plane. This structure is referred to as “orthorhombic”. We hrst use 
the PBE+U+J method to calculate the total energies of these structures at zero pressure, 
which are presented in Table [Til We consider ferromagnetic and antiferromagnetic states (the 
spin arrangement for the antiferromagnetic state for each structure is explicitly shown as or¬ 
ange arrows in Fig. [TOj) . The zero of energy is chosen as the energy of the antiferromagnetic 
state of the 2H structure, which has the lowest total energy among all the structures consid- 
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TABLE 11: Total energies E with the zero point chosen as the total energy of 2H structure with 
antiferromagnetic ordering, AE = E{AF) — E[F) the energy difference between ferromagnetic 
ordering and antiferromagnetic ordering for each structure, and equilibrium volume per Mn 
atom, calculated using the PBE+f7+J method with [7 = 5 eV and J = 0.7 eV. 


structure 

2H 

AH 

cubic 

tetragonal 

orthorhombic 

ordering 

EM 

AFM 

EM 

AEM 

EM 

AFM 

EM 

AFM 

EM 

AFM 

E (meV/Mn) 

103 

0 

186 

109 

837 

775 

724 

746 

727 

754 

AE (meV/Mn) 

-103 

-77 

-62 

22 

27 

volume H (A^/Mn) 

71.6 

71.0 

66.4 

66.1 

62.9 

61.9 

65.9 

64.1 

65.8 

63.4 


ered. We find that at ambient pressnre all the perovskite strnctures of BaMnOs are abont 
0.7 eV/Mn higher in energy that the 2H strncture. For each strncture, we also calculate the 
ferromagnetic-antiferromagnetic energy difference Ai? = E{AF) — E{F). We see that with¬ 
out Mn off-center displacement {2H, AH and cubic structures), the antiferromagnetic state 
has lower energy than the ferromagnetic state. With Mn off-center displacements (tetrag¬ 
onal and orthorhombic structures), the ferromagnetic state has lower energy. This further 
indicates that Mn off-center displacements help stabilize ferromagnetism. We also calculate 
the volume per Mn atom of each structure. We notice that antiferromagnetic state of the 
cubic structure has the smallest volume 61.9 per Mn, while the antiferromagnetic 2H 
structure that has the lowest total energy has a much larger volume 71.0 A^. Face-sharing 
oxygen octahedra are in fact more closely packed than corner-sharing oxygen octahedra. 
However, in the 2H structure, between the “column” of face-sharing oxygen octahedra there 
is hollow space, which increases the volume. 

The volume comparison suggests that external pressure can stabilize the cubic per¬ 
ovskite structure over the hexagonal 2H structure as found experimentally 26|]. To test 
this hypothesis, we calculate the total enthalpy as well as total energy of the antifer¬ 
romagnetic 2H structure and antiferromagnetic cubic structure. The enthalpy difference 
AiL = H{2H) — iL(cubic) and the energy difference AE = E{2H) — £'(cubic) are shown in 
Fig. fTTK. We hnd that at a critical pressure of 20 GPa, there is a transition from the 2H 
structure to the cubic structure. Compared to the energy difference AE, we can clearly see 
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FIG. 11: Pressure dependence of the 2H structure and the cubic structure. For the 2H structures, 
we study the antiferromagnetic state; for the cubic structure, we study the G-type antiferromagnetic 
state. The results are calculated using the PBE+t/+J method with [7 = 5 eV and J = 0.7 eV. A) 
Enthalpy difference between the 2H structure and the cubic structure /S.H = H[2H) — 77(cubic) 
(red curves). Energy difference between the 2H structure and the cubic structure AE = E{2H) — 
E(cubic) (blue curves). The green line highlights the structural transition point. B) Volume 11 per 
Mn atom of the 2H structure (red curves) and the cubic structure (blue curves). 


that it is the volume difference AV = V{2H) — l/(cubic) that stabilizes the cubic structure 
via pressure. A 20 GPa pressure is experimentally achievable and has been applied on cubic 


SrMnO, 


To make the discussion complete, we present in Fig. [TTl B the volume per 


Mn of the 2H structure and of the cubic structure as a function of pressure. We see that 
because the 2H structure is more hollow than the cubic structure, it is more compressible 
(the volume decreases faster with pressure). However, up to the critical pressure, the vol¬ 
ume of the 2H structure is always larger than that of the cubic structure, which results 
in the pressure-driven structural transition. We need a few comments here: i) we verified 
that with the application of pressure, there is no magnetic transition for the 2H structure 
and the cubic structure, i.e. the antiferromagnetic ordering always has lower energy than 
the ferromagnetic ordering; ii) at the critical pressure, the Mn off-center displacements are 
suppressed, i.e. the tetragonal and orthorhombic perovskite structures are reduced to the 
cubic structure; hi) we repeated all the calculations using the sPBEsol and sPBE methods 
and hud all the results are qualitatively consistent with a critical pressure of 21 GPa for 
sPBEsol and 25 GPa for sPBE. 
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VI. CONCLUSIONS 


We have systematically studied the phase diagram of perovskite Sri-a^BaajMnOs (0 < x < 
1) as a function of Ba concentration x and epitaxial strain, as well as the phase diagram of 
hexagonal polymorphs of BaMnOs as a function of external pressure. 

For perovskite Sri_ 3 ;Ba 3 ;Mn 03 , we hnd that both increasing the Ba concentration and 
imposing epitaxial strain tend to stabilize a ferromagnetic-ferroelectric state. In the vicin¬ 
ity of the phase boundary between a ferromagnetic state with a larger polarization and an 
antiferromagnetic state with a smaller polarization, either applying a magnetic held or an 
electric held can induce a hrst-order phase transition in polarization or magnetization, and 
thus realize an ehective giant magneto-electric coupling. However, the details of the phase 
diagram, in particular the critical Ba doping and critical epitaxial strain strongly depend on 
the choice of exchange correlation functional. We hnd that among the three havors of PBE 
exchange correlation functionals (sPBEsol, PBE-I-17-I-J and sPBE-|-17efr), sPBEsol least fa¬ 
vors the ferromagnetism and predicts the largest critical Ba concentration and largest critical 
strain, while sPBE-|-f4fr most favors the ferromagnetism and therefore predicts the smallest 
critical Ba concentration and smallest critical epitaxial strain (above some Ba concentration, 
epitaxial strain is not even needed to stabilize the ferromagnetic-ferroelectric state). We also 
show that a high pressure of over 20 GPa can stabilize the perovskite structure of BaMnOs 
over its hexagonal polymorphs. 

Perovskite Sri_^Baa;Mn 03 with x > 0.5 has not yet been synthesized in bulk or grown 
in thin him form [7|. We hope that our theoretical predictions on multiferroic properties of 
perovskite Sri_ 3 ;Baa;Mn 03 could stimulate further experiments. 
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